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Abstract. The nature of the strong absorption features in the 
white dwarf GD 229 has been a real mystery ever since it was 
found to be magnetic in 1974. All attempts to explain the spec- 
trum by line components of hydrogen failed. With the first sets 
of newly calculated line data for He I in a strong magnetic field 
we could identify most of the absorption structures in the spec- 
trum with stationary line components in a range of magnetic 
fields between 300 and 700 MG. This is much lower than previ- 
ously speculated and is comparable to the highest fields found 
in hydrogen rich magnetic white dwarfs. The reason for the 
large number of spectral features in GD 229 is the extreme ac- 
cumulation of stationary components of transitions in He I in a 
narrow interval of field strengths. 
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1. Introduction 

Magnetic fields from about 40 kG to 1 GG have been detected 
in about about 50 (2%) of the 2100 known white dwarfs (Mc- 
Cook & Sion 1998). A list of all currently known magnetic 
white dwarfs is found in Jordan (1997). Magnetic fields in the 
low field range (<20MG) in both hydrogen and helium rich 
magnetic white dwarfs can be detected by a relatively sim- 
ple line pattern. Hydrogen line components have been identi- 
fied in many objects, but until recently, helium was detected 
unambiguously only in the magnetic white dwarf Feige7. Us- 
ing He I line data calculated by Kemic (1974), Achilleos et al. 
(1992) found that the spectrum could be well reproduced with 
a dipole of strength 35 MG displaced by 0. 15 white dwarf radii 
from the stellar center. Since the Kemic data were calculated 
by regarding the magnetic field as a small perturbation to the 
Coulomb interaction and are only valid up to about 20 MG, 
the Kemic tables had to be extrapolated for the Feige 7 anal- 
ysis. The only other definite helium rich objects with fields of 
about 15-20 MG were discovered by Reimers et al. (1998) in 
the course of the Hamburg ESO survey. 

Send offprint requests to: S. Jordan, jordan@astrophysik.uni-kiel.de 



Magnetic white dwarfs with field strengths >100MG con- 
taining hydrogen could be interpreted ever since the numeri- 
cal calculations of energy level shifts and transition probabil- 
ities for bound-bound transitions by groups in Tubingen and 
Baton Rouge (Forster et al. 1984; Rosner et al. 1984; Henry 
and O'Connell 1984, 1985). With the help of these data Green- 
stein 1984 and Greenstein et al. 1985 could identify the hitherto 
unidentified features in Grw +70°8247 with stationary compo- 
nents of hydrogen in fields between about 150 and 500 MG. 

Since the magnetic field on the surface of a white dwarf 
is not homogeneous but e.g. better described by a magnetic 
dipole, the variation of the field strength from the pole to the 
equator (a factor of two for a pure dipole field) smears out most 
of the absorption lines at larger magnetic field strengths. How- 
ever, a few of the line components become stationary, i.e. their 
wavelengths go through maxima or minima as functions of the 
magnetic field strength. These stationary components are visi- 
ble in the spectra of magnetic white dwarfs despite a consider- 
able variation of the field strengths. 

About 80 % of all known white dwarfs have nearly pure hy- 
drogen atmospheres (spectral type DA). However, since most 
of the remaining stars have helium rich atmospheres, we would 
also expect a significant number of magnetic white dwarfs to 
belong to the spectral type DB in which He I lines are ob- 
served. Therefore it could be expected that the few magnetic 
white dwarfs with unidentified spectral features, that cannot 
be explained by the line data for hydrogen, are helium rich. 
The most famous example is GD 229, where Swedlund et al. 
(1974), Greenstein et al. (1974), Landstreet & Angel (1974), 
Liebert (1976), Greenstein & Boksenberg (1978), and Schmidt 
et al. (1996) found strong absorption features in the optical, 
infrared, and UV. Angel (1979) already proposed that the ab- 
sorption bands in this star are due to stationary components of 
hydrogen or helium. 

The basic difficulty in calculating energy levels and os- 
cillator strengths for arbitrary field strengths lies in the fact 
that the Coulomb potential has a spherical symmetry whereas 
the magnetic field induces a cylindrical symmetry. This pre- 
vents a separation of variables and together with the nonlinear 
character of the interactions makes even the numerical solu- 
tion a difficult problem. In particular the intermediate regime 
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(0.01 < 7 = 5/2.35 • 10 9 G < 1), in which the majority 
of magnetic white dwarfs is found, is very complicated, since 
neither the magnetic nor the Coulomb interaction dominates 
and we therefore encounter a so-called nonperturbative regime 
where none of the interaction terms can be treated by pertur- 
bation theory. In the case of the two electron system He I the 
situation is even more difficult, since the number of degrees of 
freedom is enhanced significantly and the electron-electron in- 
teraction introduces a third competitive force which makes He I 
a complex system with a rich variability of the spectrum with 
changing field strength. The particular challenge of a compari- 
son of atomic data with observational ones lies in the fact that 
the energies of many excited states have to be known with a 
high relative accuracy (usually < 10~ 4 ) and for a large number 
of field strengths in order to identify the stationary transitions. 

Until this year the status of electronic structure calculations 
on He I was such (see Braun et al. 1998, Jones et al. 1997, Ruder 
et al.1994 and references therein) that a conclusive comparison 
was impossible and the occurence of He I in the atmosphere of 
GD 229 had to be considered as a pure speculation as it was the 
case with respect to hydrogen in Grw +70°8247 until 1984. 

Therefore, several alternative explanations have been pro- 
posed. Engelhardt & Bues (1995) have tried to explain the reg- 
ular almost periodical structure of the GD229 spectrum by 
quasi-Landau resonances of hydrogen in a magnetic field of 
2.5 GG. Ostreicher et al. (1987) speculated that two of the fea- 
tures could be due to intersections of hydrogen components in 
a field of about 25-60 MG. 

Now, the calculations performed in Heidelberg (Becken & 
Schmelcher 1998) have closed that gap and for the first time 
precise data for a large number of energy states have become 
available. In this paper we will show that most of the absorp- 
tion features in the optical and UV spectrum of GD 229 can be 
identified as stationary line transitions of He I. 

2. Helium in strong magnetic field 

The fixed-nucleus electronic Hamiltonian of the helium atom 
in a strong magnetic field oriented along the z— axis reads as 
follows 

n = £ i=ll2 (^P? + &BL Xi + ^B 2 {x 2 + yf) 

+ ^ BS *>) + T7T-7TI (1) 

where the sum includes the one particle terms in the order: 
field-free kinetic, orbital Zeeman, diamagnetic, Coulomb at- 
traction and spin Zeeman energies. The last term represents 
the repulsive electron-electron interaction. Effects of the fi- 
nite nuclear mass can easily be taken into account by us- 
ing the corresponding scaling relations (Becken & Schmelcher 
1998). According to its symmetries the Hamiltonian (1) pos- 
sesses four independent conserved quantities: the total spin 
S 2 , its z— component S z , the z— component L z of the to- 
tal spatial electronic angular momentum and the total spatial 
z— parity H z (parity is a combined symmetry of the previous 
ones) which will be used to label the electronic states in the 



form n 2S+1 M^~^ z where n specifies the degree of excita- 
tion within a given symmetry subspace. Calculations are per- 
formed separately for each subspace with given quantum num- 
bers. Since the spin part of the eigenfunctions of the Hamilto- 
nian (1) are the usual spin singlet and triplet functions we ad- 
dress in the following exclusively the method of investigation 
for the spatial part of the wave functions. 

The key ingredient of our method for the ab initio calcu- 
lation of the electronic structure of atoms in strong magnetic 
fields is a basis set of one particle functions which takes on the 
following appearance 

$i(p,(j),z) = p n ">z n *>exp(-a i p 2 - /3 iZ 2 ) exp (im^) (2) 

where n Pi = \rrii\ + 2k t , n Zi = w Zi +2li with ki,k = 0, 1,2, ... 
{cti , pi } are nonlinear variational parameters; Tr Zi is the z parity 
of the one particle function. The basis set (2) has its precursor 
in a more general one (Schmelcher & Cederbaum 1988) which 
has successfully been applied for the investigation of the elec- 
tronic structure of small molecules in strong magnetic fields 
(Kappes & Schmelcher 1996). 

Due to the parameters {at, fa} our basis set is extremely 
flexible which represents a major advantage in the presence 
of an external magnetic field. This allows us to accurately de- 
scribe the changes in the electronic wave functions with chang- 
ing field strengths. The values of the variational parameters 
{(Xi(B) , 0i(B)} are determined by a nonlinear optimization 
procedure which is performed for the spectrum of both hydro- 
gen and He II in a magnetic field. Typically 150 one particle 
functions of type (2) are optimized in order to describe 15 hy- 
drogenic states with a high accuracy. The optimization proce- 
dure has to be accomplished for each field strength separately 
(see Becken & Schmelcher 1998 and references therein). 

The basic idea for solving the time-independent 
Schrodinger equation belonging to the Hamiltonian (1) is 
then to construct out of the above optimized one particle 
functions a set of spatial two electron configurations {tpi} 
being symmetric (spin singlet) or antisymmetric (spin triplet) 
with respect to the interchange of the electronic coordinates 
and respecting the above-discussed symmetries. We represent 
the exact wave functions by a linear combination of two 
electron configurations tyj = 2~)i Cji^Mri; r 2)- By virtue of 
the variational principle we arrive at a generalized eigenvalue 
problem whose eigenvalues and eigenvectors are excellent 
approximations to the exact eigenenergies and wave functions. 
We remark that all matrix elements of the Hamiltonian (1) with 
respect to the basis afunctions (2) can be evaluated analytically 
resulting in a complicated series of higher transcendental func- 
tions whose numerically stable and efficient implementation 
required a number of different techniques. 

The typical number of two-electron configurations {ipi}, 
i.e. the dimension of the full Hamiltonian matrix to be diago- 
nalized, is of the order of 2000-4000 depending on the symme- 
tries. On a powerful workstation 10 - 20 h CPU are necessary to 
build up and diagonalize the Hamiltonian matrix for each field 
strength and symmetry subspace. Extensive calculations have 
been performed in the whole range < 7 = 5/2.35 • 10 9 G < 
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100, for 30 electronic states the energy could be determined 
to an accuracy of at least 10~ 4 . For vanishing magnetic quan- 
tum number Al = many excited singlet and triplet electronic 
states with positive and negative z— parity have been deter- 
mined. For nonzero angular momentum M = +1, — 1 both sin- 
glet and triplet states have been studied for positive z— parity. 
This gives us insight into an essential part of the spectrum of 
the He I atom with varying field strength and allows us in par- 
ticular to decide whether He I can be present in the atmosphere 
of magnetic white dwarfs through an extraction of the station- 
ary components of the transitions in the corresponding regime 
of field strengths. 

3. Interpretation of the spectrum 

The strongest evidence that GD 229 does not contain hydrogen 
came from a HST spectrum taken by Schmidt et al. (1996): 
There is no evidence for components of Lyman a which do 
not vary very much over a very large range of magnetic field 
strengths. The lsrj — > 2p_i transition has a maximum at 
1341 A and should always be visible even at extremely high 
field if significant amount of hydrogen were present. 

For the first time it is now possible to compare the spec- 
tral features of GD229 directly with the newly calculated 
wavelength tables for the line components of He I in a strong 
magnetic field (Becken& Schmelcher 1998). Since oscillator 
strengths are not yet available we could not directly compare 
the spectrum with theoretical models (see e.g. Jordan 1992). 
Therefore we compared the positions of stationary line com- 
ponents with the observed absorption bands. In total about 150 
line components have been calculated of which about 40 show 
stationary behaviour at different field strengths. 

If one looks at the helium line data up to 7 = 100 (cor- 
responding to 235 GG) it is rather striking to see that the 
majority (26) of all stationary components lies in the range 
300<-B/MG<700 ( see Fig.[l]; note that the wavelengths were 
determined by spline interpolation in a still relatively crude 
grid). So far only 10 grid points exist for the range of magnetic 
field strengths plotted in Fig. [j] 

In practically all cases the minimum wavelength of the 
large number of features seen in the optical, IR, and UV spec- 
trum of GD 229 can be attributed to the wavelength minima of 
the helium components (listed in Table [l]). This explains why 
the spectrum of GD 229 has such an extremely large number of 
absorption features. 

Moreover, there are no clear contradictions between the 
predicted positions of absorptions and the observations. The 
only shortcoming of the current analysis is the fact that the 
strongest features at 4000 - 4200 A and at about 5280 A can 
only be attributed to two line transitions (2 1 + — » 2 1 _ , 
2 X + — > 2 1 (— 1) + ). Since the identification at other wave- 
lengths is convincing we believe that additional stationary 
line components are hidden in the subspaces of bound-bound 
transitions that have not been calculated yet (four subspaces 
1 (±1) _ , 3 (±1)~, and transitions with magnetic quantumnum- 
bers \m\ > 1). 
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Fig. 1. The spectrum of GD 229 (Schmidt et al. 1996) is com- 
pared with the positions of He I components which are station- 
ary in the plotted range of field strengths. The stationary wave- 
lengths are marked and labeled with a number corresponding 
to the identification in Tab. [l| Note that additional stationary 
line components which are not yet theoretically calculated may 
contribute 



In principle the wavelength maxima (also shown in the fig- 
ure) of the 2 X 0+ -> ^(-l)- and 2 J 0+ -> 2 1 (-1)~ tran- 
sitions would fit well to the position of the two strongest ab- 
sorption. However, the corresponding field strengths of about 
50 MG would be much too low compared to the result found 
for the other features. Such a large spread of magnetic field 
strengths is also rather improbable since the observed features 
are relatively sharp. 

An alternative interpretation would be that hydrogen is ad- 
ditionally present in the spectrum. The strong feature at 4124 A 
in the spectrum of Grw +70°8247 (which has a polar field 
strength of 320 MG, e.g. Jordan 1992) caused by the hydro- 
gen transition 2sq — > 4/o could also contribute to the strongest 
absorption feature in GD 229. However, in this case we would 
expect several other absorption troughs, e.g. at 5830 A (2s — > 
3po) which are not present in GD 229. This is also compati- 
ble with the finding by Schmidt et al. (1996) that no Lyman a 
components are present in the UV spectrum. 
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Table 1. Components of He I which are stationary at 300 < 
B/MG < 700 and additionally two maxima at « 50 MG. Note 
that the wavelengths (and corresponding field strenghts) of the 
minima and maxima are interpolated in a relatively crude grid. 
The error range in A was estimated from interpolations of dif- 
ferent orders 



~~No. Component Zero field trans. (A/ A) |(B/MG) 



1 


2 1 0+ 


— > 2 1 0~ 


2 1 5 


— ► 3 1 Po 


4294 ±8 


338 min 


2 


2*0+ 


-> 3 1 0- 


2 1 So 




3260 ± 6 


425 min 


3 


2*0+ 


-> 4 X 0- 


2 1 S 


-4^0 


2943 ±10 


455 min 


4 


2*0+ 


-> 5 1 0- 


2 1 5o 


-» 5^0 


2802 ± 5 


477 min 


5 


3 X 0+ 


-> 3 1 0- 


3^0 


-»4 1 F 


10692 ± 40 


379 min 


6 


3*0+ 


-» 4 1 0" 


3^0 


-^Po 


7927 ± 40 


406 min 


7 


3*0+ 


-•• 5 1 0- 


3^0 


-> 5^0 


6984 ± 40 


415 min 


g 


1 3 0+ 


-> 2 3 0~ 


2 3 S 


^3 3 P 


j ru j _i_ \j 


535 min 


9 


1 3 0+ 


-> 3 3 0~ 


2 3 S 


^4 3 P 


2504 ± 3 


658 min 


10 


1 3 0+ 


^ 4 3 0- 


2 3 S 


- 4 3 P 


2295 ± 3 


689 min 


1 1 


1 3 0+ 


-> 5 3 0~ 


2 3 S 


^5 3 P 


2200 ± 5 


705 min 


12 


2 3 0+ 


-v 3 3 0~ 


3 3 5 


^4 3 P 


9350 ± 6 


573 min 


13 


2 3 0+ 


^ 4 3 0- 


3 3 S 


^4 3 P 


6990 ± 10 


643 min 


14 


2 3 0+ 


-> 5 3 0~ 


3 3 S 


^5 3 P 


6184 ± 10 


674 min 


15 


2*0+ - 


*2 1 (-1) + 


2 1 S - 


-» 3 J P-i 


5285 ± 3 


| 50 max 


15 


2*0+ - 


*2 1 (-1) + 


2 1 S - 


-v 3 x P-i 


4812 ± 30 


258 min 


16 


2 X 0+ - 


*3 1 (-1) + 


2 1 S - 


-> 4 1 F_i 


4125 ± 3 


56 max 


16 


2 1 0+ ~* 3 1 (-1) + 


2 1 S'o 4 J F_i 


3417 ± 10 


429 min 


17 


2*0+ - 


>4 1 (-1) + 


2 1 5 - 


-> 4 x P_i 


3012 ± 10 


466 min 


18 


2*0+ - 


*5 1 (-1) + 


2 1 5 - 


-> 5 x P-i 


2825 ± 3 


523 min 


19 


3 X 0+ - 


*4 1 (-1) + 


3^0 - 


-> 4 x P_i 


8449 ± 150 


445 min 


20 


3*0+ - 


*5 1 (-1) + 


3^0 - 


-v 5^-1 


7160 ± 20 


499 min 


21 


1 3 0+ - 


>2 3 (-l) + 


2 3 S - 


^3 3 P-i 


3582 ± 2 


415 min 


22 


1 3 0+ - 


*3 3 (-l)+ 


2 3 S - 


-> 4 3 P_i 


2615 ±5 


629 min 


23 


1 3 0+ - 


+ 4 3 (-l)+ 


2 3 S - 


-> 4 3 P_i 


2338 ± 5 


689 min 


24 


1 3 0+ - 


> 5 3 (-l) + 


2 3 S - 


- 5 3 P_i 


2208 ± 3 


689 min 


25 


2 3 0+ - 


*4 3 (-l)+ 


3 3 S - 


-> 4 3 P_i 


7413 ± 40 


658 min 


26 


2 3 0+ - 


> 5 3 (-l) + 


S 3 ^, - 


-> 5 3 P_i 


6296 ± 40 


658 min 



4. Discussion and future prospects 

For the first time we have been able to prove that GD 229 is 
indeed a helium rich (DB) white dwarf at a rather high mag- 
netic field strength. /^From the magnetic field strengths corre- 
sponding to the extrema of the wavelengths of the line com- 
ponents we conclude that the range of magnetic fields is about 
300-700 MG, so that it would be compatible with a dipole field 
with a polar field strength between 600 and 700 MG. With our 
analysis we could reject previous ideas that the atmosphere is 
hydrogen rich. 

In the near future the rest of the line data for He I will be 
calculated. We expect that some of the missing line components 
will have wavelength extrema near the strong features at about 
4100 and 5300 A. 

A detailed modelling of the spectrum and polarization of 
GD 229 will be possible as soon as oscillator strengths become 
available; this will enable us to determine the magnetic field 
structure which is impossible by our wavelength analysis. 



With the new line data for He I we will also try to inter- 
pret other magnetic white dwarfs with unexplained features: 
e.g. HE 121 1-1707 (Reimers et al. 1994, Jordan 1997), and four 
other newly discovered objects from the Hamburg ESO Survey 
(Reimers et al. 1998). 

Another good candidate for such an analysis is LB 1 1 146B 
for which a convincing analysis of the spectrum with hydro- 
gen line components arrived at an approximate field strength 
of about 670 MG (Liebert et al. 1993, Glenn et al. 1994). Addi- 
tionally a very strong and broad feature at about 5600 A was 
present in the spectrum which cannot be due to hydrogen. 
Therefore Glenn et al. speculated that helium might be respon- 
sible. However, in our current set of line data no stationary line 
component of He I was found at about 670 MG. 
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